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We numerically investigate the propelled motions of a Janus particle in a periodically phase- 
separating binary fluid mixture. In this study, the surface of the particle tail prefers one of the binary 
fluid components and the particle head is neutral in the wettability. During the demixing period, the 
more wettable phase is selectively adsorbed to the particle tail. Growths of the adsorbed domains 
induce the hydrodynamic flow in the vicinity of the particle tail, and this asymmetric pumping flow 
drives the particle toward the particle head. During the mixing period, the particle motion almost 
ceases because the mixing primarily occurs via diffusion and the resulting hydrodynamic flow is 
negligibly small. Repeating this cycle unboundedly moves the Janus particle toward the head. The 
dependencies of the composition and the repeat frequency on the particle motion are discussed. 


I. INTRODUCTION 

Self-propelled motions of micro- and nano-particles 
have attracted much interest from a wide range of view¬ 
points. They will provide us with important applications, 
such as nanomachines and drug delivery m- Recent, 
focus has been on their collective dynamics because they 
are very fascinating in a growing field of non-equilibrium 
physics, i.e., active matter EHIS]. Self-propelled parti¬ 
cles use some energy or nutrients to generate the self¬ 
propulsion force. For example, biological molecules, such 
as ATPase and myosin, convert chemical energy to me¬ 
chanical motion through chemomechanical coupling El- 
In non-biological systems, Marangoni effect can induce 
spontaneous motions of liquid droplets [mus]. 

Janus particles, which have heterogeneous surface 
properties, are often employed as artificial self-propelled 
systems [ni ng. For example, the self-propelled mo¬ 
tions are modelled by asymmetric nanoparticles partially 
coated with platinum. The catalytic decomposition of 
hydrogen peroxide, which occurs selectively on the Pt- 
surface, drives the nanoparticles mu 01111 HO]- Interfa¬ 
cial phoretic effects [mill] are another possible mecha¬ 
nism of micro-swimmers m- Jiang et al. demonstrated 
that a Janus particle can create an asymmetric temper¬ 
ature gradient around it in a defocused laser beam. The 
induced gradients lead to spontaneous drift motions of 
the Janus particles [7]. Also, local heating by illumina¬ 
tion light induces active motions of Janus particles in 
a binary mixture of lower critical solution temperature 

OEU. 

Besides the thermo- and diffusio-phoresis motions [H 
EH El]; Janus particles can move in phase-separating 
binary mixtures because of the coupling between the 
wetting and phase separation [HJ ES] • The phase sep¬ 
aration of binary fluid mixtures has been well studied 
EZHUj. During the later stage of the phase separation, 
the domain patterns grow with time. The particles in 
the phase-separating mixtures are trapped in one of the 
phases or at the interfaces. Even when Brownian motions 
and external forces are absent, the particles move with 


the resulting coarsening of the domain patterns pnUHHj - 
Because the Janus particles have asymmetric wettability, 
we expect that their motions also become asymmetric. 
The direction of the asymmetric motion will more or less 
depend on the particle direction. Janus particles with 
two distinct wettabilities are occasionally used as surfac¬ 
tants to stabilize the phase-separated domains [Mh36) . 
In the final stage of the phase separation, the particle 
motions will be frozen. 

In this article, we demonstrate a possible mechanism 
of spontaneous motions of a Janus particle in period¬ 
ically phase-separating mixtures. By continually vary¬ 
ing the temperature or pressure slightly above and below 
the transition point, one can cause periodic processes of 
phase separation and mixing |571 - HU| . By resetting the bi¬ 
nary mixtures to the one-phase state, we expected that 
we could continuously propel the particle. Here, we ex¬ 
amine this expectation by means of numerical simula¬ 
tions. The dependencies of the particle motion on the 
average composition and the duration of the cycle are 
discussed. 

In Sec. jn] we explain our numerical model, which is 
based on fluid particle dynamics method [m EH . Nu¬ 
merical results are shown and are discussed in Sec. Imj 
We summarize our study and discuss some remarks in 

Sec.lWj 


II. MODEL AND SIMULATION 
A. Free energy functional 

We consider the case, in which a spherical Janus par¬ 
ticle is suspended in a binary mixture. The surface of 
the particle is heterogeneous in wettability. For numeri¬ 
cal simulations, dealing with the particle in a continuous 
manner is convenient. We express it with a smooth shape 
function as mnni: 

^ {tanh — ^ + l| ■ (1) 
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Here, r is the coordinate in a lattice space and R is 
the position of the particle in an off-lattice space, a 
is the radius of the particle and d represents the width 
of the smooth interface. In the limit of d —?► 0, '0 is 
unity and zero in the interior and exterior of the particle, 
respectively. We also define the surface distribution as 

A = 

The free energy function comprises two parts as |321 

113111], 
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( 2 ) 


The first part J^mix is the mixing free energy given by 
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where is the local concentration of the first compo¬ 
nent of the mixture. T is the temperature and we set the 
Boltzmann constant to unity, vq is the molecular volume. 
The coefficient of the gradient term C is related to the 
interface tension and is of the order of HD- 
is the Bragg-Williams type of the mixing free energy as 


/bw('/') = (l)\n(j)+ {1- (f)) ln(l -(/))+ X(^(l - (f)), (4) 


where x is the interaction parameter between two fluid 
components. Under the mean field approximation, x = 2 
and (j) = 0.5 give the critical point. The third term in the 
integrand of Eq. ^ is introduced to prevent the solvent 
from penetrating into the particle [H] . Xp and are the 
control parameters. 

Jsur is the surface free energy, which is given by 


FIG. 1. (a) A sketch of a spherical particle with a heteroge¬ 

neous surface. The orientation is described by a unit vector 
n. We set Wo = — 1 and Wi = 1 in Eq. ([^, so that we have 
W = 0 and W = — 2 at the two poles. We refer to the pole 
with lU = 0 as “head” and the other as “tail”. The particle 
head is neutral in wettability and the tail prefers component 
A. (b) A typical phase diagram of a binary fluid mixture. We 
change the x parameter by using a square wave function of 
t. Phase separation occurs during the period with X = Xd- 
However, when x = Xm, the phase-separated domains are 
mixed. 


Here p is the material density. In this work, we assume 
that all the materials have the same density. S is the 
viscous stress tensor, which is given by 

^ = T]{Wv + {Wvf}. ( 8 ) 

In the spirit of fluid particle dynamics (FPD), we assume 
that the viscosity p depends on the particle distribution 

as [mill], 

r]{r)=po + Ar]iljRir), (9) 
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= —[ drijjs{r){(j>{r) - (j)s}W{r - R,n), (5) 
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where n is the unit vector toward the orientation of the 
particle. W represents the heterogeneity of the particle 
wettability. As illustrated in Fig. Ba) , we use 

W{r - R,n) = Wo + Win- (6) 

\r-R\ 


where rjo is the viscosity of the binary fluid and Arj is the 
viscosity difference between the solvent and the particles. 
In the limit of Ap —)■ oo, the particles will behave as solid 
particles [12] ■ p is a part of pressure, which imposes the 
incompressibility condition: V • = 0. 

f is the force field stemming from the particle interac¬ 
tions and is given by, 


tljR{r) dF 
n dR’ 


( 10 ) 


where Wq and Wi are the material constants for the wet¬ 
ting. If W <0, the component of larger (p tends to wet 
the surface [ 23 ] 11 ]' la the early stages of phase separa¬ 
tion, the surface of lU 7 ^ 0 largely influences the pattern 
formation (see below). 


B. Time development equations 

The hydrodynamic equation for the flow field v is given 
by 

p(4 + „.v). 

= + -^ + ^ ^ 9 “ W +VS. ( 7 ) 


where U is the effective volume defined as fl = 
J dripjf^{r). It is approximated as U « 47ra^/3 in three- 
dimensional systems, g originates from the torque acting 
on the particle, and is given by 
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The particle motions are caused by the hydrodynamic 
flow and their kinetics are described as 




V = ^ y driPji{r)v, 

dr'tjj^{r)W x x n. 


( 12 ) 

(13) 

(14) 
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Furthermore, the time development equation of the con¬ 
centration field is 

+ +V-C, (15) 

where -F(V') is the kinetic coefficient, and in which we 
set = Lo(l ~ to eliminate the flux inside the 

particles. Lq is the kinetic coefficient of the bulk mix¬ 
ture. C represents the thermal fluctuation satisfying the 
fluctuation-dissipation relation. In this model, the dif¬ 
fusion flux does not contribute to the particle motion 
directly. In each state, the total free energy, including 
the kinetic energy, should decrease with time. Its tem¬ 
poral change is described in Appendix]^ Our model can 
be applied to many-particle systems, where Janus parti¬ 
cles behave as surfactants We hope that we will 

report our studies on them elsewhere in the near future. 


C. Numerical simulations 

We numerically solve the above equations using the 
Maker and Cell method with staggered grid [33] . We dis¬ 
cretize the space by d and set vq = and C = d'^. Also, 
we set the particle radius to a = 6d. The simulation box 
is a three-dimensional system (64^) with periodic bound¬ 
ary conditions. The time increment is 0.005 toj where 
to is a typical diffusion time defined by to = d^T/LQ. 
For the wettability, we set Wq = — 1 and Wi = 1 in 
this study. Hence, the head of the particle is neutral 
in the wettability and the tail prefers the component of 
large cj) (see Fig. [^a)). Hereafter, we describe the more 
and less wettable components as A and B, respectively. 
The viscosity parameters in Eq. ([^ are rjo = pL^/T and 
A 77 = 49770 . In Eq. (1^, we set ^o to be equal to the av¬ 
erage composition, and Xp = 20. Because Reynold 
numbers in colloidal systems are very small, we iterate 
to integrate Eq. Q without updating i?, n, and (j) until 
\p{d/dt + V ■ V)iq becomes less than I0“^?7oTo/(3od^)- 
The intensity of the thermal fluctuation is given by 
(C^(r,^)Co(^'^tO) = 0-05TL{tp)6{r — r')S{t — t')5ij, where 
i and j stand for x, y and z. 


D. Periodic phase separation 

To induce periodic phase separation, we change the x 
parameter uniformly in space using square wave functions 
of time (see Fig. ^b)). and to denote the durations 
for mixing and demixing, respectively. We also define 
tt = tm + td- For simplicity, because we consider mix¬ 
tures near the phase-separation point, we assume that 
the other physical parameters are constant, independent 
of X- In durations of mtt < t < mtt + tmi X is fixed 
at Xm below the critical point. Here, m is an integer. 
The phase-separated domains will be mixed during these 
times. In mtt + tm < t < (w + I)tt, we retain x at 


X = Xd, above the critical point, so that phase sepa¬ 
ration proceeds in the bulk. Throughout the simula¬ 
tions presented in this article, we fix the x parameters 
to Xd = 2.7 and Xm = 1-3. In the demixing periods 
with Xd, the coexistence concentrations of A- and R-rich 
phases are (j) — 0.893 and 0.107, respectively. The corre¬ 
lation length, also referred to as the interface thickness, 
is = ■\/C'/(x — Xc) — 1.20fi. However, during the mix¬ 
ing periods with Xm, the correlation length is given by 
Cm = x/C'/{2(Xc - x)} = 0.845d. 


III. RESULTS AND DISCUSSIONS 
A. Domain patterns 

Figure shows snapshots of the domain patterns in 
the case of tm = td = 100 to- The snapshots were ob¬ 
tained at t = 200 to and t = 1000 to, which correspond to 
the times when the first and fifth demixing periods have 
just finished, respectively. The average concentrations, 
{4>), of component A are (a) 0.3, (b) 0.5, and (c) 0.7. In 
the demixing periods with Xd, the volume fractions of 
the A-rich phase are (a) 24.6%, (b) 50% and (c) 75.4%, 
respectively. The light blue surfaces represent the isosur¬ 
faces of (f) = 0.5 and the dark blue sphere represents the 
Janus particle. At t = 0, the particle is directed toward 
the z-axis, he., n{t = 0) = (0,0,1). In each case, we 
observed that the particle follows the background flows, 
which are caused by the interface tensions during the 
phase separations. The mechanisms of the particle mo¬ 
tions are discussed below. 


B. Particle trajectories 

Figurej^shows typical trajectories of the Janus particle 
from t = 0 to t = 10^ to in th® binary mixtures of different 
concentrations. The time intervals are tm = td = 500 to- 
The trajectories show that the particle tends to move di¬ 
rectionally toward the particle head. Because the particle 
surface has a heterogeneous affinity to the components, 
the phase separation proceeds asymmetrically around the 
particle. This asymmetry of the phase separation dynam¬ 
ics may cause the directional motion of the Janus particle. 
Interestingly, the particle moves toward its head in all the 
mixtures. Figurej^also indicates that the trajectories are 
not completely straight. The degree of the directionality 
and the particle speed depend on parameters such as the 
average concentrations and the time intervals. 

Figures l^a) and (b) show the temporal changes in the 
trajectory length and velocity toward the head. Since 
the particle orientation, n, changes with time, the tra¬ 
jectory length d|| and the velocity V\\ toward the particle 
orientation are calculated as 


d\\(t)= fmt')dt', 

Jo 


(16) 
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FIG. 2. Snapshots of a Janus particle in periodically phase- 
separating binary fluids. The average concentrations of the 
more wettable phase are (a) (0) = 0.3, (b) {(j)) = 0.5, and (c) 
(0) = 0.7. 



FIG. 3. The trajectories of the Janus particle in the periodi¬ 
cally phase separating mixtures from t = 0 to f = 10“^ to- The 
trajectories are projected on the (a) x-z and (b) x-y planes. 
At t = 0, the particle is directed toward the 2 -axis. The time 
intervals are td = tm = 500 to. 


= (17) 

Positive and negative values of represnt the forward 
and backward motions of the particle position, respec¬ 
tively. In Fig. 1^ the time intervals are set to tm = = 

10^ to- Here, each curve was obtained from one simu¬ 
lation run. We simulated nine average concentrations 
from {(j)) = 0.1 to 0.9. In the mixtures of {(j)) <0.1 and 
{(p) > 0.8, we did not observe any drastic motion of the 
particle; hence, their curves have not been included in 


Fig.S Also, we did not plot the curves for the mix¬ 
tures of {(p) = 0.4 and {(p) = 0.6 because they essentially 
demonstrate behaviors similar to those with the symmet¬ 
ric mixture {(p) = 0.5. 



FIG. 4. The time developments of (a) the particle displace¬ 
ments d|| and the particle velocity V\\ toward the particle head. 
The time intervals are td = tm = 10^ to- The average concen¬ 
trations are {(p) = 0.2, 0.3, 0.5 and 0.7. The hatched regions 
indicate the mixing periods. 

In Fig. I^a), the trajectory lengths indicate stepwise 
motions. The particle is almost fixed in the mixing peri¬ 
ods. However, during the demixing periods, the particle 
shows forward displacements. By repeating these cyclic 
motions, the particle continuously propels in a periodi¬ 
cally phase-separating binary mixture. The displacement 
in each cycle is of the order of the particle diameter. Al¬ 
though the onsets of the motion in the demixing peri¬ 
ods are not clearly seen in Fig. the detailed analyses 
indicate that the particle does not start moving simul¬ 
taneously with the quenching into the demixing states. 
It moves most largely after a certain incubation time t,, 
which is discussed later. 

The cyclic behaviors are also clearly displayed in the 
particle velocity. After the initial incubation time in each 
demixing period, the particle velocity shows large posi¬ 
tive values. After this transient deterministic motion, the 
velocity decreases gradually with some fluctuations. In 
particular, it can have negative values in more symmetric 
mixtures, with {(p) = 0.5. A similar stepwise motion is 
observed in a system where a Janus particle with metal¬ 
lic surfaces is trapped at a liquid-air interface [20]. The 
particle’s stepwise motion is due to spontaneous cyclic 
bursts of bubbles. In our system, the stepwise motion is 
due to the controlled changes of the interaction parame¬ 
ter. 

Figures]^ a) and (b) show the temporal changes of the 
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average-concentration differences and the velocity inten¬ 
sities. They are calculated as 


{Ap^) J dr{p- {P)f, 

(18) 


(19) 


where fit is the system volume. In the early stages of the 
demixing periods, the phase separation starts via diffu¬ 
sion of the components. As shown in Fig. the hy¬ 
drodynamic flow develops simultaneously with the phase 
separation, and thus, it is small during the early stages 
[46l Wj\ . Because the particle is transported by the back¬ 
ground flow, the incubation time of the particle motion 
ti in Fig. [^corresponds to the duration of the early stage 
of the phase separation, te- Here, te depends on the aver¬ 
age concentration. This is because the growth rate of the 
concentration held depends on the average concentration. 
Inside the spinodal regime, the phase separation proceeds 
via spinodal decomposition. As the average concentra¬ 
tion approaches the spinodal points ((^) = 0.5 ± 0.255 
for Xd = 2.7), the growth rate is decreased to zero. How¬ 
ever, in the binodal regime, the phase separation occurs 
via the nucleation of the droplets of the minority phase. 
The nucleation rate is also decreased to zero as {(p) ap¬ 
proaches the equilibrium concentration. In both the pro¬ 
cesses, more symmetric mixtures are more unstable and 
the durations of the early stages are shortened. This may 
suggest that the particle moves faster in the symmetric 
mixtures because the incubation time, ti, during which 
the particle is at rest, is reduced. We define the duration 
of the early stage te as -I- 1 ^ -|- te)) = (A^g^)/ 2 . 

Here, A^eq is the concentration difference in the equilib¬ 
rium demixing state. From the simulation results, we ob¬ 
tain te = 34.6 to for {(p) = 0.5, te « 74.0 to for {p) = 0.3 
and 0.7, and te « 292 to for {p) = 0.2. These are indi¬ 
cated by the arrows in Fig. [^a). 

Figs, j^a) and (b) show plots of the averaged speeds 
toward the particle head, (V||), and perpendicular to it, 
(Vj_), as a function of the time interval. Here, we set 
fm = for simplicity. The parallel and perpendicular 
velocity in the m-th cycle are defined as 

Vllm = - (20) 

H Jmtt 

V^m = - |V(t) X nit)\dt. (21) 

J mtt 

From these, the averaged velocities (Vj|) and (Vj_) are 
obtained with averaging M = 10 cycles as 

1 “ 

(Vx) = Vxm, ( 22 ) 

m—1 

where X refers to A =|j and T. The error bars in 
Figs, j^a) and (b) represent the standard deviations of 
V||j„ and Vim- The ratio of (VI) to (V||) is plotted in 
Fig. j^c). 



FIG. 5. The time developments of (a) the average- 
concentration differences (A^^) and the velocity intensity 
(n^). The time intervals are td = tm = 500to. In (b), the 
curve for (p) = 0.2 is magnified tenfold. The average concen¬ 
trations are (p) = 0.2, 0.3, 0.5 and 0.7. The hatched regions 
indicate the mixing periods. 

In Fig. ^a), each curve of (V||) is non-monotonic with 
maxima. The maxima peaks suggest that we can choose 
efficient time intervals for propelling the particle. As dis¬ 
cussed above, the hydrodynamic flow has not developed 
yet in the early stage of the demixing periods. In the 
cases of small time intervals, the phase-separating times 
are too short for the hydrodynamic flows to develop suffi¬ 
ciently; hence, the particle is not dragged, implying that 
we have to maintain the system in the phase-separated 
state till the conclusion of the early stage te- However, 
for large time intervals, the average speed along the ori¬ 
entation becomes low. In the mixtures of (p) = 0.2, 0.3, 
and 0.7, the parallel speed becomes largest at approxi¬ 
mately td — te- Thus, larger time intervals are not needed 
for propelling the particle with high speeds. However, in 
the symmetric mixture {{p) = 0.5), the parallel speed be¬ 
comes largest around tu « 200 to, which is larger than the 
duration of an early stage tg. This is attributed to large 
background flows emerging from other domains near the 
particle, as discussed below. 

As has been shown in Fig. we observe the particle 
fluctuations, which are characterized by the perpendic¬ 
ular motion (V(l) in Fig. §;b). The particle changes its 
orientation and the resulting direction of the particle mo¬ 
tion. Thus, similar to the (Vj_) changes, the temporal 
changes of the orientation are also considered as a mea¬ 
sure of the particle fluctuation. Fig. [^a) demonstrates 
the autocorrelation of the orientation vector. It is calcu¬ 
lated as 

1 P^max 

N{t) =- / dt'n{t') • n{t + t'), (23) 

^max Jo 
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FIG. 6. Plots of the averaged velocity (a) toward the particle 
head (V^|) and (b) perpendicular to the orientation (Ki), with 
respect to the time intervals. The arrows in (a) indicate the 
durations of the early stage te- The average concentrations 
are (</)) = 0.2, 0.3, 0.5 and 0.7. The error bars represent the 
standard deviation of the particle motions, (c) The ratio of 
the perpendicular speed (Fl) to the parallel speed (Vj|). 

where we set tmax = 5 x In Fig- [^a), we plot 

the autocorrelation for (0) = 0.3 as a typical example. 
They decrease with time indicating that the memory of 
the orientation is gradually lost. Figure a) shows that 
the decay rate is increased with increasing The au¬ 
tocorrelations for different (</>) behave in the same way. 
In Fig. I^b), we plot the change rate of the orientation 
which is calculated for M = 10 as 

M-l 

= M E (24) 

m—0 

n^ = ^{n(im+l)tt)-n{Tntt)}. (25) 

As shown in Figs. I^b) and Bb), both (VI) and 
((n^))^/^ are large for large values of fd) in contrast to 
those for small In the demixing periods, the phase 
separation proceeds throughout the bulk. The fluctua¬ 
tions of the particle motions and orientation stem from 
the hydrodynamic flows accompanied by the spontaneous 
growth of the domains surrounding the particle. These 


background flows are independent of the particle; thus, 
they disturb the particle motion and change its orienta¬ 
tion. Therefore, if we use large time intervals, the particle 
motion is likely to deviate from the straight line along the 
initial orientation n{t = 0). As suggested in Figs. and 
the ratios (V1)/(V||) and ((n)^)^/^ increased slightly 
with td- Also, in this sense, the large time interval is not 
preferred for the controlled propulsions. 

Figure l^a) shows that the propelling speed is higher 
in the symmetric mixtures {{(j)) = 0.5) than in the asym¬ 
metric mixtures. However, the fluctuations of the par¬ 
ticle motion are also large as shown in Figs. I^b) and 
m- Therefore, the symmetric mixtures are not suit¬ 
able for keeping the straight line motions. The particle 
moves largely; however, its motion easily loses the di¬ 
rectionality with time. Furthermore, in the asymmetric 
mixtures, the fluctuations of the particle motion and the 
orientation are relatively small for the preferred td; thus, 
the asymmetric mixtures are more suitable to control the 
particle motions. We consider the mechanisms of the pro¬ 
pelled motions in the two types of asymmetric mixtures 
separately. 


C. The wettable component-rich mixture 

Figure shows the patterns of the evolutions of the 
concentration and flow field in the mixture of (</>) = 0.7. 
As shown in Fig. the particle moves rather straightfor¬ 
ward. In the early stage of the phase separation during 
the demixing periods, the A-rich phase wets the half¬ 
portion of the particle and a wetting layer is formed on 
it. Because of the directional fluxes of the A-component 
toward the surface, the A-component is depleted near 
the outside of this first wetting layer. This accumulated 
layer structure is similar to the oscillating profiles of the 
concentration field near a flat wall or a particle with a 
homogeneous surface [47l |48] . To compensate for the de¬ 
pleted region, non-spherical droplets of the B-rich phase 
are formed as shown in Fig. [^b). The droplets grow 
near the surface via coagulation and coalescence. Be¬ 
cause the processes of the coagulation and coalescence 
occur asymmetrically near the particle tail, the associ¬ 
ated hydrodynamic flow likely pushes the particle toward 
the particle head. This process is shown in Figs.j^c) and 
(d). During the phase separation, many small droplets 
of the B-rich phase are also formed in the bulk. They 
grow to their typical sizes with time via coalescence and 
coagulation, or evaporation and condensation processes 
[27l[29j. Around coalescing droplets, other hydrodynamic 
flows are induced, which disturb the particle motion. The 
disturbing flows are relatively weaker; thus, the particle 
motion remains straightforward, in contrast to that in 
the symmetric mixtures. 

The mixture of {(j)) = 0.8 should be phase-separated 
at equilibrium, when y = Xd(= 2.7). However, because 
the nucleation rate is very small, the time intervals we 
employed (td < 10^ to) were not enough to induce the 
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phase separation. Thus, we did not observe any motion 
of the Janus particle in the mixture of {(p) > 0.8. 



tih 



FIG. 7. (a) The autocorrelation function of the orientation 

N{t). The average concentration was constant at {p) = 0.3, 
whereas the time intervals changed, (b) The dependence of 
the change rate of the particle orientation on the 

time interval. 



Ml to = 0 Ml to = 40 



7^A;/Io=100 Atlto=200 


FIG. 8. Snapshots of typical pattern evolution around the 
particle in the A-rich mixtures {{(j>) — 0.7). The arrows show 
the flow held. The red and blue domains are the A- and B-rich 
phases, respectively. The time intervals are = tm = 300 to- 


D. The wettable component-deficient mixture 

The motions in the i?-rich mixtures are highly directed 
toward the particle orientation as shown in Figs. and 
1^ Figure shows the pattern evolutions of the concen¬ 
tration and flow field in the mixture of {(p) = 0.3. In the 
early stages of the phase separation, a cap-shaped domain 
of the A-rich phase emerges from the wettable portion of 
the particle surface. It covers a large amount of the in¬ 
terfacial area, although its volume is not so large. Then, 
it tends to change its shape to a sphere for reducing the 
interfacial energy after the early stage of the phase sep¬ 
aration. This process induces a pumping hydrodynamic 
flow around the tail of the Janus particle and the result¬ 
ing pumping flow pushes the particle toward the head as 
shown in Figs.j^c) and (d). The decay time of the local¬ 
ized pumping flow, is estimated as th ~ ''ya/c, where 
a is the interface tension. 

The mixtures of (p) = 0.3 and 0.7 have the same sta¬ 
bility for the phase separation in the bulk. Because the 
tail of the particle prefers component A, this asymmetry 
leads to the difference in the particle motions between the 
(p) = 0.3 and the (p) = 0.7 cases. Figure|^a) shows that 
the highest parallel speed in in the (p) =0.7 mixture is 
approximately twice of that in the (p) = 0.3 case. How¬ 
ever, Figure l^c) indicates that the ratio {Vx)/{V\\) in 
the {p) = O.Tcase is also approximately twice of that in 
the {p) = 0.3 case at the maximum parallel speed. Also, 
Figure |^a) shows that the standard deviations of Vj_m 
for the (^) = 0.3 mixture are smaller than those for the 
(p) = 0.7 case. The particle is suggested to move more 
steadily and smoothly in the (p) = 0.3 mixture. Thus, 
we conclude that the JJ-rich mixtures are more preferred 
to induce more straight motions. 

The stability for the bulk phase-separation in the mix¬ 
ture of (p) = 0.2 is the same as that in the (p) = 0.8 
case. However, we have not observed any motions in the 
(p) = 0.8 mixture, although the particle moves straight¬ 
forward in the {p) = 0.2 case. This difference suggests 
that the domain formation in the mixture of (p) = 0.2 is 
attributed to the heterogeneous nucleation at the parti¬ 
cle surface [13] ■ As the average concentration approaches 
the binodal line, the thermal nucleation rate in the bulk 
is strongly decreased |13|. However, the rate of heteroge¬ 
neous nucleation on the wetting surface is large enough 
to induce it during our demixing periods. Because the 
number of the surrounding droplets is decreased, parti¬ 
cle motion becomes more straightforward in more asym¬ 
metrically H-rich mixtures (see the case of (p) = 0.2 in 
Fig.|§). 

In both cases, the hydrodynamic flow around the parti¬ 
cle tail pushes the particle toward the head. In this sense, 
our particle motion may be categorized to a pusher in the 
active matter field m- However, more detailed analyses 
on the flow pattern are required before this conclusion 
can be derived. 

















E. Roles of mixing periods 


Figures |^a) and (b) indicate that the particle does 
not show any large changes of the position and orienta¬ 
tion during the mixing periods. This is because diffusion 
dominates the mixing but it does not contribute to the 
particle motion. However, here, we should note that the 
mixing periods are very important for resetting the bi¬ 
nary mixtures. As noted above, the hydrodynamic flow 
caused by the interface tension is large in the demixing 
periods, whereas it is small in the mixing periods. This 
difference in the hydrodynamic flows present between the 
mixing and demixing periods induces the continuous pro¬ 
pelled motion of the Janus particle. 

In the above simulations, we set tm = td for simplic¬ 
ity. However, considering whether this mixing interval is 
sufficient to reset the binary mixtures is important. The 
characteristic length i of the phase-separated domain in¬ 
creases with time algebraically. In the symmetric mix¬ 
ture, a bicontinuous pattern is formed and the domain 
grows obeying [25] 

^(t') « cat If]. (26) 


(a) 






A/Z/o = 0 

o 

11 ^ 
o 

<1 

l! Jis'> 



At/to =100 

At/to =200 


FIG. 9. Snapshots of typical pattern evolutions around the 
particle in the R-rich mixtures (((^) = 0.3). The arrows show 
the flow field. The red and blue domains are the A- and R-rich 
phases, respectively. The time intervals are td = tm = 300 to- 


However, in the asymmetric mixtures, the minority phase 
forms droplets. The droplets grow with time via coales¬ 
cence and coagulation as [271111] 

e{t) « c'iTtfrjy/^. (27) 

Here, c and c' are non-dimensional numbers, which de¬ 
pend on the volume fractions. Then, the characteris¬ 
tic length at the ends of the demixing periods would 
be given by £{td)- If the volume fraction is quite small, 
the droplets grow via nucleation and growth mechanism, 
and the above scaling relation (Eq. @) is replaced by 
i{t) (X {L^atY^^. Because its growth exponent is the 
same as that of Eq. , we consider only cases described 
by Eq. ^2n\ above. 

In the mixing periods, these phase-separated domains 
should be dissolved into a homogeneous state via diffu¬ 
sion. The diffusion time is estimated as tdif ~ ^(td)^/JJm, 
where = 2Lo(Xsp~Xm)/7^ is the diffusion constant at 
X = Xm and Xsp is the interaction parameter at the spin- 
odal point. When the time interval of the mixing period 
is longer than the diffusion time, the system can be reset 
for the next demixing period. However, if tdn ^ tm: the 
mixtures show non-steady states [57H1D] and the particle 
would not move directionally. The condition tdif ^ tm 
is rewritten as for the droplet 

patterns and tm ^ {a/rjYt'^/Dm for the bicontinuous 
patterns. If we set td ^ (T/ry)^ZJ“^, we have tdif ^ td 
in the asymmetric mixtures. Therefore, long annealing 
times for the mixing are not required. The total time 
interval ft = tm + td can be reduced to it — td- 


IV. SUMMARY AND REMARKS 

We proposed a possible mechanism of propelled mo¬ 
tions of Janus particles in periodically phase-separating 
binary mixtures using numerical simulations. Because 
the particle has an asymmetric surface in wettability, 
phase separation proceeds heterogeneously around the 
particle. The resulting asymmetric hydrodynamic flow 
leads to directional motions of the particle. After a usual 
one-step quench, the directional motion slows down and 
will stop eventually. By changing the interaction param¬ 
eter, the phase-separated system recovers to a one-phase 
mixing state. Therefore, under periodic changes of the 
interaction parameter, we can induce a continuous mo¬ 
tion along the particle orientation. 

We found that the particle propels more directionally 
in asymmetric binary mixtures. In symmetric mixtures, 
the hydrodynamic flow emerging from other domains sur¬ 
rounding the particle is so large that it disturbs the direc¬ 
tional motions of the Janus particle. We also found that 
the propelling speed and the directionality depend on the 
frequency of the change in the interaction parameter. In 
the cases of short time intervals, the hydrodynamic flow 
does not develop well and it cannot drive the particle. 
However, in the cases of long time intervals, the hydro- 
dynamic flows from the surrounding domains disturb the 
directional motion. We can efficiently move the particle 
in the intermediate time intervals, which are comparable 
to the characteristic time of the spinodal decomposition 
or nucleations in the early stage of the phase separation. 
The pumping hydrodynamic flow that is localized around 
the particle decays with the relaxation time th- is of 
the order of rja/a. The interval of the demixing period 
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should be larger than ti + th- However, the large time 
interval for the demixing periods leads to the loss of the 
directionality as discussed above. Thus, the most efficient 
time interval would be td ~ ti + th- The displacement of 
the particle in each cycle is of the order of the particle 
diameter. Then, the maximum particle speed could be 
increased toH~a/(ti + th)- 

Unfortunately, the simulations performed in this study 
are limited owing to the numerical costs. We need to 
deepen our understanding of the propelled motions of 
such a Janus particle, and we hope that we will present a 
report on them in the near future. We make some critical 
remarks to improve our study as follows. 

(1) In this article, we show only the simulations with 
the particle diameter 2a = 12d. Here, d is comparable to 
the correlation length; hence, our particle is rather small. 
In actual phase-separating mixtures, such small particles 
would show drastic Brownian motions and the directional 
motions we discovered might be smeared out. Our pre¬ 
liminary simulations with larger particles indicated that 
the displacement of the particle in each cycle is of the 
order of the particle size. In other words, they suggest 
that the particle speed can be increased linearly with its 
size by employing appropriate temporal changes of the 
interaction parameter. 

(2) The wettability of our particle changes smoothly 
on the surface (see Eq. (©)• However, an actual Janus 
particle usually has two distinct surfaces and the wetta¬ 
bility changes abruptly at the equator. Instead of Eq. (§, 
we performed some simulations with an alternate surface 
function given by 


Wir 


R, n) = Wq + Wi tanh 



r-R\ 

\r-R\) 


(28) 


Here, dh is introduced to avoid the singularity at the 
equator of the Janus particles. A particle described 
by Eq. (28) with small d^ would behave more realisti¬ 
cally. The preliminary simulations using Eq. (28) with 
dh = 0.05 demonstrated essentially similar results as 
those using Eq. Q. This similarity can be derived from 
the fact that our particle motions are induced by the sec¬ 
ondary effect of the phase separation, i.e., the asymmet¬ 
ric growth of the hydrodynamic flow around the particle. 
Thus, we consider that our findings are robust for the 
heterogeneous surface pattern. To improve the efficiency 
of the particle motions, simulations with a variety of the 
surface structures would be interesting. 

We consider that the propelled motions in this article 
are robust also for the details of the free energy function. 
We numerically confirmed that the particle moves in the 
same way in binary mixtures described by the Ginzburg- 
Landau free energy instead of Eq. Q. 

(3) The easiest method to realize our findings in ac¬ 
tual systems would be a direct observation of them with 
an optical microscope equipped with a temperature con¬ 
trol hot stage [40]. However, changing the temperature 
with high frequencies might be experimentally difficult 
because the thermal diffusion constant is finite. To in¬ 


duce a high speed propulsion, large differences of the tem¬ 
perature from the binodal point are preferred for both 
the mixing and demixing periods. However, for deep 
quenches in the demixing periods, the incubation times 
for the phase separation would be reduced up to micro¬ 
scopic timescales. 

The pressure control m is considered to be an alter¬ 
native method to induce the continuous propulsions. The 
combination of the temperature control and illumination- 
induced phase separation would also be able to induce 
periodic phase separation with high speeds |53J [21] . In 
this study, we employ only the square waves of the in¬ 
teraction parameter to induce periodic phase separation. 
We should perform more simulations with other types of 
wave functions to find more efficient propulsion schemes. 

(4) In our model, particle motion is caused by the hy¬ 
drodynamic flow. However, studying the behaviors of a 
Janus particle in solid mixtures where = 0 would also 
be interesting. In such mixtures, the particle motion is 
caused by the thermodynamic forces (oc \7SR/6(j)). Be¬ 
cause the coarsening behavior of the phase separation 
domain pattern depends on the fluidity of the mixtures 
(see Eqs.(|2^ and {2^), whether our findings are appli¬ 
cable to the solid mixture is not trivial. 

(5) Because the intensity of the hydrodynamic flow is 
proportional to l/rjo (see Eq. 0). we expected that the 
particle speed is also proportional to l/rjo. However, our 
simulations with different solvent viscosities indicated 
that Vii is approximately proportional to I /770 hi symmet¬ 
ric mixtures, whereas it is approximately independent of 
rjQ in asymmetric mixtures (not shown here). This is be¬ 
cause the rate-limiting process of the phase separation 
in asymmetric mixtures is the nucleation and growth of 
the minority phase; their rates are dominated by the dif¬ 
fusion constant, not the solvent viscosity . Thus, the 
results reported in this article are not quantitatively uni¬ 
versal. They will depend on the system parameters such 
as the solvent viscosity and interface tension. Although 
we consider that the propelled motion in periodic phase 
separation can be qualitatively realized in actual systems, 
we have to perform more simulations with wide ranges of 
the system parameters. 

For simplicity, we also assumed a constant viscosity tjq 
in the surrounding fluids. However, the solvent viscosity 
is generally a function of the local composition. As noted 
above, the particle speed depends on the solvent viscosity 
in symmetric mixtures; thus, studying the influences of 
the viscosity difference would also be interesting. 
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Appendix A: Time development of total free energy 

The total free energy including the kinetic energy is 
considered. Using Eq. ([^, it is given by 

Q{(l),v,R,n) = F+ \ f drpv"^. (Al) 


Its time development is 
dr 


-g = 

dt 


1 

— 0 + pU • V 

dF ■ dF . 

m' 


After some calculations, we obtain 

2 


dt 


g= dri-LlW 


,dF 

Jcj) 


(A2) 


— Vu : H + f ■ V 


• (V X u) + ( p + -pv'^ ) V • u 


+ 


(dF ■ OF . 

ISr " 


(A3) 
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where we neglect the contribution of the system bound¬ 
ary. For the force / and torque g fields, we obtained the 
following equations in our scheme. 


J drf ■ V 


1 


(9J- 


-- / dri){r - R)— ■ v{r 


V 


dR 


1 OR /■ , ,, ^ 

~VdR' - R)v{r) 


dR 

dR 


R, 


(A4) 


^ J drij^ij — H){V X u(t')} X n 
Then, we finally obtain 

= +5’)(v» + (v«)t|, 

(AC) 


dn 


J drg • (V X u) 

'ip{r — R) 


= — I dr 


V 


n X 


dR 

dn 


{V X v{r)} 


where we assumed V • u = 0. Thus, we confirmed that Q 
decreases with time only via diffusion and viscous dissi¬ 
pation. 



